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A comprehensive study of the electronic properties of monoclinic M0O2 from both an experimental 
and a theoretical point of view is presented. We focus on the investigation of the Fermi body and the 
band structure using angle resolved photoemission spectroscopy, de Haas- van Alphen measurements, 
and electronic structure calculations. For the latter, the new full-potential augmented spherical wave 
(ASW) method has been applied. Very good agreement between the experimental and theoretical 
results is found. In particular, all Fermi surface sheets are correctly identified by all three approaches. 
Previous controversies concerning additional hole-like surfaces centered around the Z- and B-point 
could be resolved; these surfaces were an artefact of the atomic-sphere approximation used in the 
old calculations. Our results underline the importance of electronic structure calculations for the 
understanding of M0O2 and the neighbouring rutile-type early transition-metal dioxides. This 
includes the low-temperature insulating phases of VO2 and Nb02, which have crystal structures 
very similar to that of molybdenum dioxide and display the well-known prominent metal-insulator 
transitions. 
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I. INTRODUCTION 

The metal-insulator transitions of the early transition- 
metal oxides have been attracting a lot of attention for 
decades. In particular, since the first report by Morin, 1 
much interest has centered about the cornerstone ma- 
terials V2O3 and VO2. This is mainly due to the fact, 
that these compounds display first order phase transi- 
tions with very narrow hystereses of only few Kelvin and 
strong changes of the conductivity of several orders of 
magnitude, making them promising candidates for appli- 
cations. Interestingly, the transitions of both compounds 
are accompanied by characteristic structural changes, 
which initiated a discussion about the driving forces. Dis- 
pute centered about the question whether the transitions 
are driven predominantly by structural instabilities of the 
parent corundum and rutile structure, respectively, or 
else by strong electronic correlations. Nowadays, the im- 
portance of the latter is widely accepted. However, de- 
spite intense work, the questions concerning the origins of 
the respective transitions have not yet been satisfactorily 
answered. 

In VO2, the metal- insulator transition occurs at 340 K 



and is connected with a structural change from the ru- 
tile structure to a distorted monoclinic structure. The 
latter is characterized by (i) pairing of the metal atoms 
within chains parallel to the rutile c-axis and (ii) their 
lateral antiferroelectric zigzag- like displacement Elec- 
tronic states near the Fermi energy are mainly of V 3c? 
t2 g character. They fall into the d\\ band, which me- 
diates V-V overlap along the metal chains, and the re- 
maining, bands4 At the transition, splitting of the 
dy band and upshift of the bands due to increased 
metal-oxygen overlap produce a finite band gap. Contro- 
versial discussions addressed the origin of the dy splitting, 
which was assigned either to the metal-metal dimeriza- 
tion or to increased electronic correlations resulting from 
the reduced screening; by the el electrons State 
of the art electronic structure calculations gave strong 
hints at a structural instability but were not able to 
reproduce the insulating gap due to the shortcomings 
of the local density approximation (LDA) . 10 ! 11 ' 12 Only 
recently, LDA+DMFT calculations have demonstrated 
that the metal-insulator transition may be regarded as a 
correlation-assisted Peierls-like transition^ 

Remarkably, the structural distortions occurring at 
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the transition of VO2 are characteristic of most of the 
neighbouring d , d 2 , and d 3 transition-metal dioxides. 
In particular, the metallic oxides M0O2, WO2; TCO2, 
and a-Rc02 all display the above-mentioned metal-metal 
pairing as well as the lateral displacement and crystal- 
lize in the same monoclinic structure as low-temperature 
VO2 M- This crystal structure is displayed in Fig.[Tl which 
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FIG. 1: (Color online) Monoclinic crystal structure of M0O2. 
The dimerized metal pairs are highlighted by thick bonds. 

allows to identify both the pairing and the lateral dis- 
placement of the metal atoms. 

In contrast to the above-mentioned metallic d 2 and d 3 
members, the d 1 compound Nb02, like VO2, undergoes 
a metal-insulator transition (at 1081 K) and a simultane- 
ous structural transition from rutile to a distorted vari- 
ant having a body-centered tetragonal (bet) lattice. 15 ' 16 
However, despite the differences in long-range order, the 
local deviations from the rutile structure are the same as 
in VO2, i.e. the niobium atoms dimcrizc and experience 
lateral displacements at the transition. 17 ' 18 ' 19 

Despite of the strong interest in the metal-insulator 
transition of VO2, attempts at understanding the origin 
of the common structural characteristics are very rare. 
Indeed, any theory aiming at the metal-insulator transi- 
tion of VO2 would be incomplete without proper account 
of the local structural distortions of all the rutilc-related 
transition-metal dioxides. In order to close the apparent 
gap, we have started a systematic investigation especially 
on VO2, NbC"2, and M0O2 already a decade ago, which 
by now provides a unified picture. 11 ' 12 ' 20 ' 21 

Interest in M0O2 was especially motivated by the 
metallic nature of this material, the reduced localization 
of the Ad electrons as compared to the 3d electrons of 
VO2, and the rather small effective mass, which is simi- 
lar to the free electron mass . 22 ' 23 ' 24 These facts point at a 



rather weak influence of electronic correlations. Indeed, 
electronic structure calculations for M0O2 show a large 
bonding-antibonding splitting of the dy bands i 11 ' 20 Since 
the occupied dy bands lie far below the Fermi level, a 
possible influence of electronic correlations in these bands 
would not lead to an enhanced effective mass or suscep- 
tibility. As a consequence, M0O2 can be regarded as 
a "model system" for studying the origin of the struc- 
tural characteristics of the early transition-metal diox- 
ides. In addition, an extensive examination of the elec- 
tronic structure of M0O2 could help clarifying open ques- 
tions of the above-mentioned structurally related sys- 
tems, in particular of the d 1 members NbC>2 and VO2. 

In a previous work, we have presented a compre- 
hensive theoretical study of the electronic properties of 
MoC^. 11 ' 20 Electronic structure calculations using the 
ASW method, which is based on density functional the- 
ory and the local density approximation (LDA), were 
used to analyze the electronic states involved in low- 
energy excitations. In addition, emphasis was put on 
investigating the origin of the small but distinct devia- 
tions of the monoclinic crystal structure from the par- 
ent rutile structure, which for this material mainly con- 
sist of a metal-metal dimerization parallel to the rutile 
c-axis. In short, we were able to show that the mono- 
clinic structure results from what we called an embedded 
Peierls-like instability. This instability affects the quasi 
one-dimensional Mo 4dy bands, which mediate the metal- 
metal overlap along chains parallel to the rutile c-axis 
and which are embedded in a background of the remain- 
ing three-dimensionally dispersing Mo 4d e^, electrons, 
which are also designated as the ir* electrons. The valid- 
ity of this scenario had been demonstrated by comparing 
the electronic properties calculated for the monoclinic 
structure with those of a properly symmetrized artificial 
rutile structure. In addition, we found nearly perfect 
agreement of our calculated results with the outcome of 
ultraviolet photocmission spectroscopy and x-ray absorp- 
tion spectroscopy apart from an underestimation of the 
dy band splitting by about 1 eV due to the well known 
shortcomings of the LDA. 11 ! 20 

This scenario has also been shown to explain the metal- 
insulator transitions of the d 1 members VO2 and NbC>2, 
which likewise are subject to Peierls-like instabilities of 
the quasi one-dimensional V 3dy bands in the embedding 
backgrounds of the ir* electrons! 12 ' 21 As in M0O2, these 
instabilities cause distortions of the rutile structure by 
the characteristic metal-metal dimerization. Due to the 
particular position of the Fermi energy in the d 1 com- 
pounds, the Pcicrls-likc instabilities are then accompa- 
nied by the observed metal-insulator transitions. At the 
same time, the opening of the optical band gaps is sup- 
ported by the lateral antifcrroclcctric displacements of 
the metal atoms, which cause an upshift of the tt* states. 

Continuing our study of M0O2, we concentrate in the 
present work especially on the Fermi surface and the band 
dispersions. For the first time, a comprehensive investi- 
gation of the full Fermi body of M0O2, as determined 
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both experimentally and theoretically, will be presented. 
In particular, we turn to a comparison of the theoretical 
results with angle-resolved photoemission spectroscopy 
(ARPES) and de Haas- van Alphen data (dHvA). In do- 
ing so, we apply a new implementation of the augmented 
spherical wave method, which goes beyond the version 
used in our previous work by taking into account the full 
potential rather than the approximate muffin-tin poten- 
tial. The resulting changes of the electronic states are 
small but important and lead to a much improved agree- 
ment of the theoretical and experimental data. 

While published photoemission studies on M0O2 fo- 
cused on the angle-integrated analysis of the occupied 
bands ; 25 ' 26 ' 27 ' 28 a systematic investigation of the band 
dispersions and the Fermi body by angle-resolved pho- 
toemission spectroscopy seems to be still missing. In 
addition, existing magneto-transport measurements ad- 
dressed only fractions of the Fermi surface . 22 i 23 ' 24 ' 29 ' 30 ' 31 
It is the purpose of the present work to make up for the 
apparent lack in experimental data and to present a com- 
prehensive study of the electronic structure of Mo02- 

An important difference between dHvA and ARPES is 
that the former requires very low temperatures whereas 
the latter does not. This is crucial for the early 
transition-metal dioxides, some of which have a metal- 
insulator transition and thus do not show a Fermi sur- 
face at low temperatures. In the present case we have 
the opportunity to compare ARPES and dHvA because 
Mo0 2 is metallic at low T. It turns out that the ARPES 
and dHvA results agree well, implying that ARPES will 
be useful to study the Fermi surfaces of the d 1 members 
V0 2 and Nb0 2 as well. 

The paper is organized as follows: While Sec. [Til gives 
an overview on the computational method and the the- 
oretical results, Sec. IIIII turns to a description of the 
ARPES and dHvA measurements and their results. The 
theoretical and experimental findings are subject of a 
comparative discussion given in Sec. [TV] and Sec. [Vj sum- 
marizes our results. 



II. ELECTRONIC STRUCTURE 
CALCULATIONS 

A. Computational Method 

The calculations are based on density-functional the- 
ory and the local density approximation. 32 ' 33 They were 
performed using the scalar-relativistic implementation of 
the augmented spherical wave method (see Refs. [HI! 
and references therein). In the ASW method, the wave 
function is expanded in atom-centered augmented spher- 
ical waves, which are Hankel functions and numerical so- 
lutions of Schrodinger's equation, respectively, outside 
and inside the so-called augmentation spheres. In order 
to optimize the basis set, additional augmented spherical 
waves were placed at carefully selected interstitial sites. 
The choice of these sites as well as the augmentation 



radii were automatically determined using the sphere- 
geometry optimization algorithm^ 6 - Self consistency was 
achieved by a highly efficient algorithm for convergence 
acceleration^ The Brillouin zone integrations were per- 
formed using the linear tetrahedron method with up to 
518 k-points within the irreducible wedge. 35 ' 38 

In the present work, we used a new full-potential ver- 
sion of the ASW method, which was implemented only 
very recently^ In this version, the electron density and 
related quantities arc given by a spherical harmonics ex- 
pansion inside non-overlapping muffin-tin spheres. In the 
remaining interstitial region, a representation in terms of 
atom-centered Hankel functions is used4°- However, in 
contrast to previous related implementations, it is suf- 
ficient to work without having to employ a so-called 
multiple-fv basis set, permitting a very high computa- 
tional speed of the resulting scheme. 



B. Results 

The electronic properties of M0O2 in the monoclinic as 
well as in the above-mentioned artificial rutile structure 
have already been discussed in much detail in our pre- 
vious work. 11 ' 20 For this reason, we will in the following 
concentrate on the changes coming with the considera- 
tion of the full potential rather than its muffin-tin approx- 
imation. In particular, we will investigate the changes of 
the electronic states near the Fermi energy. 

Furthermore, since a very detailed description of the 
monoclinic crystal structure and its relation to the rutile 
structure has been given in our previous work, we will not 
enter a renewed discussion here. However, we refer the 
reader to our previous publication especially concerning 
the definition of the local coordinate systems to be used 
for the analysis of the electronic states. 11 ' 20 

While a first account of the monoclinic crystal struc- 
ture was given by Magneli and Andersson^i, our cal- 
culations are based on the refined data by Brandt and 
Skapski. 42 Hence, in contrast to our previous work ; 11 ' 20 
we will exclusively deal with the measured monoclinic 
structure and not discuss our calculations for the related 
but artificial rutile structure. In addition, since there are 
no indications for the formation of magnetic moments 
in M0O2, spin degeneracy was enforced in all calcula- 
tions. The resulting partial densities of states (DOS) 
are shown in Fig. [2j They look very similar to those 
presented previously . 11 ' 20 In particular, we recognize the 
6.5 eV wide O 2p bands well below the Fermi energy Sp- 
in contrast, the Mo Ad t2 g and e g manifolds are found 
around and well above E-p. Worth mentioning is the pro- 
nounced dip of the tig partial DOS at Ep. It results from 
the Peierls-like instability of the d x 2_ y 2 states, which usu- 
ally are designated as the d\\ states and which mediate 
the metal-metal overlap along the Mo atom chains paral- 
lel to the rutile c-axis. The splitting of these bands into 
bonding and antibonding branches as resulting from this 
Peierls-like instability is clearly visible in Fig. [3J which 
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FIG. 2: (Color online) Partial densities of states (DOS). 
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FIG. 3: (Color online) Partial Mo 4d t 2s DOS. Selection of 
orbitals is relative to the local rotated reference frame, see 
Ref . [OHO- 




FIG. 4: (Color online) First Brillouin zone of the simple mon- 
oclinic lattice. 
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FIG. 5: (Color online) Electronic bands. 



displays the three partial t 2 g DOS. In representing these 
partial DOS we have used the same local coordinate sys- 
tem as in our previous work with the local £-axis parallel 
to the rutile c-axis and the local z-axes pointing alter- 
nately along the [110] and [110] direction,!^ In Fig. 
[3l we observe the strong splitting of the d x 2_ y 2 states, 
which do no longer contribute to the metallic conduc- 
tivity. Furthermore, due to the splitting, the low lying, 
bonding branch is found at the lower edge of the ti g group 
of bands and nearly separated from the higher lying d xz 
and d yz bands. The latter are commonly designated as 
the 7r* bands. 

The electronic bands along selected lines within the 
first Brillouin zone, Fig. HJ are displayed in Fig. [5j With 
the interpretation of the partial DOS at hand, we eas- 
ily recognize the O 2p as well as the Mo Ad t2 g and e g 
states in the energy intervals from —9 to — 2.5cV, —2 to 
+3eV, and +3 to +7eV, respectively To be more de- 



tailed, we find a two bands between —2 and — 1 eV, which 
is separated from all other bands. This split-off doublet is 
more easily identified in Fig. [6l where we show the near- 
Ep bands on an expanded energy scale. According to the 
previous analysis, the split-off doublet is just the bonding 
dy band. The one-dimensional behaviour of the dy bands 
is contrasted by the rather isotropic dispersion of the it* 
bands. As already pointed out in our previous work, hy- 
bridization between both types of bands is rather small. 
In contrast coupling of these states is only via charge 
conservation. For this reason, we proposed to interpret 
the splitting of the former in the course of the structural 
distortions as an embedded Peierls-like instabilit y 11 ! 20 

Concentrating especially on the Fermi body of mon- 
oclinic M0O2, we complement Fig. [5] by Fig. \7\ which 
displays cuts through the Fermi body in the ac-planc, 
the ab-plane, and perpendicular to the a-axis as well as 
a three-dimensional view of the Fermi surface sheets. 

The Fermi body is defined by the Fermi level crossings 
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FIG. 6: (Color online) Near-_BF electronic bands on an ex- 
panded energy scale. 



III. EXPERIMENTS 
A. Sample preparation and characterization 

M0O2 single crystals were grown by chemical transport 
using TeCU as transport agent M- The crystals exhibit 
specular surfaces of a size up to 3 x 3 mm 2 . The sur- 
face orientation was determined by x-ray diffraction. For 
the ARPES measurements the as-grown surfaces were 
cleaned in situ by heating to 750° C for some minutes. 
The quality and cleanness of the heated single-crystal 
surfaces were controlled by low energy electron diffraction 
(LEED) and by scanning tunnelling microscopy (STM). 
For the dHvA measurements several rods of 3 mm length 
and 1.5 mm diameter were prepared. The b-axis was the 
axis of rotation of the rods. The residual resistance ratio 
was 97 for j |j a and 144 for j || b. 



of the d xz and d yz bands. In Figs. [6] and [3 we identify 
three different Fermi surface sheets: two hole-like sur- 
faces near the Y-point (sheet 1 and 2), which occupy 9 
and 2%, respectively, of the volume of the 1st Brillouin 
zone, as well as a peanut-shaped electron-like surface cen- 
tered around the T-point and occupying 11 % of the Bril- 
louin zone volume (sheet 3). The latter sheet has two 
loopholes near the B-points, whereas it is closed along 
the direction T-A. 

As already mentioned in the introduction, use of the 
exact single-particle potential, i.e. of the so-called full po- 
tential, instead of its muffin-tin or atomic sphere approxi- 
mation leads to important improvements. These become 
obvious from comparing Fig. [5] of the present work to 
Fig. 6 of Ref. H3- Apart from a broadening of the O 2p 
bands we find distinct differences especially at the Z- and 
the B-point. To be specific, those bands, which form the 
highest occupied states at these points in the present cal- 
culation, were found just above Ep in the old calculation. 
They generated small hole-like pockets of the Fermi sur- 
face centered around these points, which were at variance 
with the experimental findings (see below). In contrast, 
in the new calculation these bands are found below Ep 
and the hole-like pockets disappear. 

Another important difference can be observed at the 
T-point. The maximum binding energy of the highest oc- 
cupied band of the new calculation (0.30 eV) agrees much 
better with experimental value of 0.21 eV (see below) 
than the binding energy of the old calculation (0.43 eV). 
Thus, by using the full-potential ASW method, the agree- 
ment between theory and experiment is considerably im- 
proved. 



B. Angle resolved photoemission spectroscopy 
(ARPES) 

Experimentally, the electronic structure E(k) was in- 
vestigated by ARPES measurements using the geometry 
and emission angles defined in Fig.[5J Due to energy con- 
servation the binding energy Eb of the photoelectron is 
determined by 

\E B \ = hv - E kin - <f>, (1) 

where hv is the photon energy and <fi the work function of 
the crystal. For a detailed description of the photoemis- 
sion process we refer the reader to the literature i 45 i 46 ' 47 i 48 
On passing through the crystal surface, the component 
of the electron wave vector parallel to the surface ky re- 
mains unchanged with 

|kii| = \^2mE kin sin<& . (2) 

The perpendicular component kj^ is, strictly spoken, not 
conserved in the photoemission process but it is usually 
possible to reconstruct its approximate value by including 
information about the final states of the photoemission 
process. If one assumes a free electron-like final state 
dispersion citenielsen03 , k± = |kj_| can be calculated as 



k± = -^2m{E km cos 2 $+\E \+4>) . (3) 

Here, Eq is the inner potential of the sample as referred 
to E-p. While the work function <f> can be determined ex- 
perimentally, the inner potential Eq is usually estimated 
from an analysis of the measured band structure. As a 
consequence of Eq. ((3]), the set of all ky vectors at a fixed 
photon energy describes a spherical surface in k-space. 

A Fermi vector kp is defined by the Fermi level cross- 
ing of a valence band in the corresponding direction of 
k-spacc. As ARPES permits a complete mapping of 
a solid's three-dimensional band structure, the position 



6 




b) 





Sheet 3 (Electrons) Sheet 2 (Y-Holes) 



Sheet 1 (Y-Holes) 





C ) «Y ( C « '• C ( «Y ) C 




Sheet 2 (Y-Holes) Sheet 1 (Y-Holes) 



Sheet 3 (Electrons) 








z z 



y ( r. . r 



z z 



r« 



y r 



Sheet 3 (Electrons) Sheet 2 (Y-Holes) Sheet 1 (Y-Holes) 



d) 




FIG. 7: Cuts through the Brillouin zone (left column) and the Fermi body a) in the ac-plane, b) in the ab-plane, and c) 
perpendicular to a as well as d) 3-dimensional view of the three Fermi surface sheets (generated using XCrysDen^). 



of the Fermi level crossings can be identified. A three- 
dimensional ARPES measurement can be performed by 
varying the emission angles and hv , causing a variation 
of k|| and k_L, respectively. 

The ARPES measurements were done at the beam- 
line SGM3 at the synchrotron ASTRfD in Aarhus, Den- 
mark. The beamline provides light in the energy range 
12 cV < hv < 130 cV with an energy resolution of 



E/AE = 150003 The UHV chamber is equipped with 
an ARUPS10 electron analyzer (Vacuum Generators) 
with an angular resolution of 0.7° and an energy reso- 
lution of 25 mcV. The analyzer is fixed on a computer 
controlled two-axis goniometer with an angular accuracy 
up to 0.2°. All ARPES data were taken at a temperature 
of 30 K. For the description of the measurement geome- 
try, instead of the polar emission angle its projections 
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FIG. 8: Measurement geometry and definition of the emission 
angles in an ARPES experiment. Here, K is the wave vector 
in vacuum. 

on the xz- and yz-planes are used (see Fig. EJ. The re- 
sulting emission angles are called $ and 9 . The direction 
of normal emission was determined by a careful optical 
alignment of the highly reflective surface using a laser. 

C. ARPES Results 

1. (100) surface 

ARPES spectra taken on the (100) surface for varia- 
tion of the emission angles <& and are shown in Fig. 

and b, respectively. At hv = 19 eV, binding energies 
up to 0.6 eV have been measured. In both graphs a sin- 
gle peak is observed, which displays a nearly parabolic 
dispersion. At normal emission ($ = and = 0), 
it reaches a maximum binding energy of 0.22 eV. With 
increasing emission angles, both $ and O, the peak ap- 
proaches the Fermi level and and crosses it. As will be 
shown later by the symmetry of the Fermi surface, the 
T-point is reached at hv = 19 eV. Hence, the maximum 
binding energy of 0.22 eV corresponds to the binding en- 
ergy at the T-point. We can thus compare the occupied 
bandwidth determined here directly with the binding en- 
ergy of the highest occupied LDA band at the T-point 
which is 0.3 eV. 

Analyzing the nearly parabolic dispersion of the 
ARPES band, effective masses of 0.876m e and 0.821m e 
can be found for variation of the emissions angles <E> 
and 0, respectively. This is in excellent agreement with 
Shubnikov-de Haas measurements of Ref. [24], which 
provided for comparable geometry an effective mass of 
0.83m e . 

2. Determination of the Fermi body 

The Fermi-level crossing of the highest occupied band, 
defining the corresponding Fermi vector, can be identi- 
fied by a strong increase of photocmission intensity at 
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FIG. 9: ARPES measurements on the (100) surface under 
variation of the emission angle a) $ and b) (hv = 19 eV). 

E-p. Indeed, ARPES can also provide a direct image of 
the Fermi surface when the photoemission intensity in 
a narrow energy window around the Fermi level is dis- 
played as a function of emission angle or photon energy 
(see Ref. [49| and references therein) . 

For a quantitative analysis of the ARPES spectra a 
constant background as determined from the energy re- 
gion of unoccupied states has been subtracted. To correct 
intensity variations caused by different emission angles, 
the spectra were normalized to this background. The in- 
tensity dependence on hv is considered by normalization 
with the average total intensity of all spectra taken at one 
certain photon energy. The work function has been deter- 
mined by VLEED measurements to cj) = 5.15 ± 0.3eVi£i 
The final electron states are supposed to be free-electron 
like. The two following methods are used to extract the 
Fermi vectors from the ARPES data: 

1. The photocmission intensity is integrated over a 
120 meV wide energy region near E-p. The maxima 
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of this integrated intensity I(k) define the Fermi 
vectors. 

2. Since the Fermi surface is determined by strong 
variations of the intensity distribution I(k), the 
gradient |Vi<J(k)| could be likewise used to de- 
fine the Fermi vectors . 52 ! 53 However, this gradient 
method generally produces two Fermi contour lines. 
Whereas the contour line on the unoccupied side of 
the Fermi surface coincides with the real Fermi sur- 
face, the other line is an artefact of this technique^ 2 - 

Both methods can be found in the literature. However, 
the gradient method often results in a better agreement 
with theoretical predictions^ 2 - 



3. Fermi surface scans perpendicular to the direction a* 

In the following, ARPES kn scans on the (100) surface 
will be presented, which allow for a Fermi surface map- 
ping of planes perpendicular to the direction a*. Such 
Fermi surface maps can be performed by variation of 
both emission angles at a fixed photon energy. They 
provide spherical cuts through k-space. Nevertheless, for 
representation reasons, flat projections of the spherical 
paths are shown. The theoretical data are produced in 
complete analogy: spherical paths in k-space are calcu- 
lated, but flat projections are shown. In order to simplify 
the discussion, we define a Cartesian coordinate system 
within the plane of the projection by k x and k yi which 
are the components of the vector ky parallel and perpen- 
dicular to the direction T-Z, respectively. This situation 
is sketched in Fig. [TOb . 




FIG. 10: (Color online) k-space cuts corresponding to kii mea- 
surements in the a) (100) surface and b) (101) surface. The 
spherical curvature of the area in ARPES measurements is 
neglected. 

In Fig. [TTh . an ARPES k|| scan through the Fermi 
body of M0O2 is displayed. The measurement has been 
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FIG. 11: (Color online) ky mapping of the a) Photoemission 
intensity 1(h) at Ep, and b) gradient of the photoemission in- 
tensity |Vi c J r (k)| at Ef, both resulting from an ARPES mea- 
surement on the (100) surface taken at a photon energy of 
hv = 19 eV. Solid (red) lines point to structures of maximal 
intensity. Points 1 to 4 mark the Fermi vectors in direction 
of k x and k v . 



performed at a photon energy of hv — 19 cV. The pho- 
toemission intensity at E-p is plotted as a function of k x 
and k y using the intensity method (method 1). Two dif- 
ferent Fermi surface structures, characterized by an en- 
hanced photoemission intensity, can be identified, which 
are marked by solid (red) lines. Most prominent is the 
ellipsoidal structure, which is located around the ori- 
gin and characterized by the points 1 to 3 with com- 
ponents k x pi = 0.175A -1 , k y F2 = — 0.175A -1 , and 
k y F3 = 0.150 A" 1 , respectively. In addition to this 
electron-like Fermi surface a much less pronounced hole- 
like Fermi surface contour is observed, albeit for nega- 
tive values of k y only. It comprises point 4 with com- 
ponent k y Fi = — 0.375 A -1 and separates electron-like 
states with k y > k y Fi from hole-like states in the outer 
region k y < k y p^. 

Figure II lb shows the photoemission intensity gradi- 
ent |Vk/(k)| as a function of k x and k y (method 2), 
based on the data shown in Fig. [TTk . In Fig. [TTb . 
the above-mentioned double structure of the gradient 
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method is nicely visible for the electron-like Fermi sur- 
face centered around the origin. The real Fermi surface 
contours, i.e. those not caused by the artefact of the gra- 
dient method, are marked by solid (red) lines. The ex- 
tent of this Fermi surface amounts to k X Fi = 0.225 A -1 , 
k y F2 = — 0.213 A -1 , and k y F3 = 0.188 A -1 , respectively. 
The hole-like Fermi surface contour for negative values 
of k y is also observed in Fig. [TTb . However, the photoc- 
mission gradient in this region is very small and, hence, 
the expected double Fermi surface structure can not be 
detected. Only the real Fermi surface contour is visi- 
ble. The extent of the hole-like Fermi surface amounts 
to kypi = — 0.438 A -1 . In general, the error of both 
methods of Fermi surface determination is defined by the 
k-resolution of the ARPES experiment and amounts to 
|Ak| w 0.025A- 1 . 
In Fig. [T^l the experimental ARPES Fermi surface con- 
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FIG. 12: (Color online) Cut through the Fermi body perpen- 
dicular to a*: ARPES and LDA Fermi surface contours. 

tours as resulting from both the use of the intensity at E-p 
and the gradient of the intensity are compared to the con- 
tours resulting from the LDA calculations. As mentioned 
above, the LDA calculations are performed on spherical 
paths, so as to account for the spherical curvature of the 
experimental plane. However, in the figures, the flat pro- 
jection on the plane perpendicular to a* is shown. Only 
the k-space region accessible by the ARPES apparatus 
is included. 

According to Fig. IT2J the LDA calculations predict an 
electron-like, elliptic Fermi surface around the r-point. 
Note that, as will be shown later, k± reaches the r-point 
right at hv m 19 cV. As a consequence, the ky mappings 
shown in Figs. [TTb and b as well as in Fig. fTJ] cover a 
plane perpendicular to a* , which includes the r-point 
(see Fig. [TUh). The theoretical findings are confirmed 
by the ARPES measurements. However, as compared to 
the LDA results, the Fermi surface contours determined 
by the intensity method are shifted towards the occupied 
states. Hence, the ARPES Fermi surface is smaller than 



the calculated one with the deviations of the measured 
Fermi vectors ranging from —5 to —32%. Here and in 
the following all deviations are given with respect to the 
calculated values. In contrast, the gradient method pro- 
duces smaller shifts in both directions, away from as well 
as towards the occupied states. In this case, the devia- 
tions of the measured from the calculated values range 
from -15 to 20%. 

Perpendicular to the T-Z direction, cuts through two 
hole-like Fermi surface sheets centered around the Y- 
point are predicted by theory (see also Fig. [7h). How- 
ever, owing to the spherical curvature of the measure- 
ment plane, observation of both hole-like Fermi surface 
contours is expected only for positive values of k y . In 
contrast, for k y < only one of these contours is cut by 
the measurement plane, which is indeed observed in the 
ARPES experiments. In the region of negative k y its in- 
ner surface can be seen in Fig. fTJl while both the outer 
part in k y < and the two hole-like structures for k y > 
are not detected. Compared to the LDA calculations, 
also the hole-like ARPES Fermi surface contour shows a 
slight shift, which, for the results obtained by the inten- 
sity method (method 1), is directed towards the occupied 
states. The corresponding ARPES Fermi surface is thus 
bigger than the calculated one, and the experimental de- 
viations from the theoretical values amount to 13 to 24 %. 
For the gradient method the ARPES contours generally 
agree well, but significant shifts with deviations between 
—2% and 13% are observed near some fractions of the 
Fermi surface. No systematic trend for the direction of 
these shifts can be discerned. 

The experimental and theoretical values of the Fermi 
vectors k x p and k y p in the plane perpendicular to a* 
through the r-point as well as the deviations of the theo- 
retical values from the experimental ones as arising from 
the gradient method are summarized in Tab. [TJ 
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I(k) 


0.175 


-0.175 


0.150 


-0.375 


\V k I(k)\ 


0.225 


-0.213 


0.188 


-0.438 


LDA 


0.207 


-0.237 


0.222 


-0.429 


Deviation 


9% 


-10% 


-15% 


-2% 



TABLE I: Magnitudes of Fermi vectors in inverse A. The ex- 
perimental Fermi vectors are taken from ky scans on the (100) 
surface (Jw — 19 eV) and analyzed using both the intensity 
maximum and its gradient. The last lines gives the respective 
deviations of the experimental result as based on the gradi- 
ent method from those arising from the LDA calculations. 
The plane perpendicular to a* through the r-point has been 
analyzed. The points 1 to 4 are marked in Fig. 1121 
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jf.. Fermi surface scans in the a* c* -plane 

A variation of kj_ can be achieved by varying the pho- 
ton energy (see Eq. ©). In the so-called angle-energy 
scans one of the emission angles as well as the photon 
energy hv are varied while the respective other emission 
angle is held constant. From this, Fermi surface scans as 
a function of kj_ and one of the parallel wave vector com- 
ponents (k x or k y ) can be obtained. In particular, angle- 
energy scans on the (100) surface under variation of hv 
and the emission angle (variation of k y ) allow a Fermi 
surface mapping of the a*c*-plane, which is spanned by 
the vectors a* and c*. With k x = ($ = 0), a cut 
through the T-point is achieved. In this situation, kj^ 
is parallel to the direction T-B, whereas k y is orientated 
perpendicular to both T-B and T-Z. The thus described 
cut through the Brillouin zone is sketched in the left col- 
umn of Fig. [7^. 

The resulting ARPES Fermi surface is given in Fig. 
[T3"k . which shows the intensity /(k) as a function of k y 
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FIG. 13: (Color online) Angle-energy scan on the (100)- 
surface: In dependence of k y and k± a) the photoemission in- 
tensity at Ef, and b) the gradient of the photoemission 
intensity at Ef, |Vk/(k)|, is plotted. Artefacts are marked by 
arrows. (Light grey (red) points: electron-like Fermi surface 
contours, dark grey (blue) points: hole-like contours) 

and k±. Both, electron- like and hole-like Fermi surface 
contours can be found, which are marked by light grey 
(red) and dark grey (blue) points, respectively. 

The inner potential, which is necessary for the cal- 



culation of k±, can be estimated by a comparison of 
the ARPES Fermi surface contours to the LDA results. 
The LDA calculations give a peanut-shaped electron-like 
Fermi surface, which is symmetric with respect to the 
T-point (see Fig. [7^,). This shape is reflected by the 
ARPES measurements, which likewise display the waist 
of a peanut of the electron-like Fermi surface. As already 
mentioned above, the center of this structure, which is 
just the T-point, is reached at hv m 19 cV. Using the 
parabolic free-electron like final state dispersion (third 
Brillouin zone), an inner potential of Eq f=a — 7eV is 
found. This value is comparable to the inner potential 
defined by the bottom of the valence band in LDA cal- 
culations and ARPES measurements {Eq^oa ~ — 9cV). 

Fig. [T3"b displays the gradient of the intensity, 
|Vic-T(k)|, again as a function of k y and k±, based on 
the data shown in Fig. [T3a . The typical twofold struc- 
ture of real Fermi surface contours and the technical arte- 
facts can be observed. The real Fermi surface contours 
are marked by points. Despite the normalization of the 
ARPES data, there remain some variations of the pho- 
toemission intensity, which are independent of the Fermi 
surface. In Fig. 113b . those artificial structures are marked 
by arrows. They will be neglected in the subsequent dis- 
cussion. 



5. Fermi surface scans perpendicular to the direction 
(a*-c*) 

Fermi surface scans have been also performed on the 
(101) surface. If again the spherical curvature of the 
measurement plane is neglected, the corresponding k|| 
scans on this crystallographic orientation produce cuts 
through the Fermi body perpendicular to the direction 
(a* — c*) (indicated in Fig. [TUb). In particular, while 
k x is parallel to the direction T-Z, k y is perpendicular 
to both the directions T-Z and T-A. This situation is 
sketched in Fig. [TUb . 

As a matter of fact, ARPES measurements on the 
MoO2(101) surface lead to less pronounced PES peaks 
than the measurements on the (100) surface. This might 
be due to a higher roughness of the surface or an en- 
hanced density of defects. 

ARPES ky scans on the MoO2(101) surface taken at 
hv = 19 cV are displayed in Fig. [TH where the intensity 
I(k) is plotted as a function of k x and k y . Near the origin, 
two separated structures of enhanced photoemission in- 
tensity can be observed, which in the direction k y extent 
to k y Fi = —0.08 A -1 and k y p2 = 0.48 A -1 , respectively. 
In contrast, in the direction k x no Fermi vector can be 
defined. The band crossing Ep in point 1 is occupied 
by electrons in the region k y > k y p\, the band respon- 
sible for point 2 is occupied by electrons in the region 
k y < k y p2- Unfortunately, the photoemission peak creat- 
ing these Fermi surface contours can not be observed over 
the whole range of k y pi < k y < k y p2- As a consequence, 
it is not clear, whether both crossing points are generated 
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FIG. 14: (Color online) k||-mapping at the (101) surface at 
hv — 19 eV. The photoemission intensity at Ef, I(k) is plot- 
ted in dependence of k x and k y . The solid (red) lines point 
up structures of maximal intensity. Points 1 and 2 mark the 
Fermi vectors in direction of k x and k v . 
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FIG. 15: (Color online) a*c*-plane of the Brillouin zone and 
of the Fermi body: Angle-energy-scans on the (100) and (101) 
surfaces produce cuts perpendicular to the particular kx . The 
end points of the corresponding Fermi vectors are marked by 
light grey (red) and dark grey (blue) points, respectively. 



by a single band or by two different bands. If they are 
due to a single band, the region k y pi < k y < k y p 2 must 
be electron-like. Then both Fermi vectors would be part 
of the calculated peanut shaped, electron-like Fermi sur- 
face around the T-point (see Fig. \7§. In contrast, if the 
crossing points were created by two different bands, only 
one of the Fermi vectors could be part of the electron-like 
Fermi surface. 

Angle-energy scans on the (101) surface allow the in- 
vestigation of the a*c*-plane (see Fig. [7^,, left column), 
where kj_ is in the direction of Y-A and k y is orientated 
perpendicular to both the directions Y-A and Y-Z. The 
corresponding results will be included in the summary 
plot for the Fermi surface contours in the a*c*-plane to 
be presented in the following section. 



6. Fermi surface in the a* c* -plane: ARPES vs. LDA 

So far, different ARPES mappings of the Fermi sur- 
face contours of M0O2 in the a*c*-plane have been de- 
scribed. Furthermore, detailed LDA calculations for this 
plane have been done. In the following, all these data 
will be combined into a single representation. To this 
end, the ARPES results on the MoO2(100) surface and 
the theoretical findings can be combined without further 
treatment. Yet, inclusion of the measurements on the 
(101) surface is more complicated. In a combined plot 
of results from the (100) and (10-1) surface orientations, 
the respective Y points must coincide. Furthermore, since 
the normal vectors of the (100) and (101) surface enclose 
an angle of 60.4°, the perpendicular wave components 
kx of the two surfaces also have to enclose this angle. 
The situation is sketched in Fig. [15] For illustration the 
LDA Fermi surface and some exemplary cuts through 
the Brillouin zone obtained by anglc-cnergy-scans on the 
two crystallographic surfaces have been included in this 



figure. However, for each crystal surface there is an un- 
certainty of the orientation of ±3°. Accordingly, for the 
angle enclosed by the two surfaces an error of ±6° is sup- 
posed. 

Of course, in the ideal case, the ARPES findings on the 
two different crystallographic surfaces of M0O2 should be 
consistent among each other and they should be consis- 
tent with the LDA calculations. Within the experimen- 
tal errors, an adjustment of the inner potential Eq, the 
work function <fi and the angle between the perpendicular 
wave vectors is possible. The best agreement between all 
experimental and theoretical findings is achieved for an 
inner potential of Eq = — 8.5cV and a work function of 
<j> = 5 eV for the analysis of the data of the (101) surface 
and an angle of 55.4°. Using these values, we combine all 
experimental and theoretical Fermi surface results on the 
a* c* -plane into Fig. [TBI Obviously, the overall agreement 
between the experimental and theoretical findings is very 
good. Large parts of the theoretically predicted, peanut- 
shaped, electron-like Fermi surface around the T-point 
are observed by the experiment. In particular, the waist 
of the "peanut" is correctly confirmed by all ARPES 
measurements. Moreover, the two hole-like Fermi sur- 
face contours around the Y-point can also be observed in 
the ARPES measurements. Especially at k y f=a —0.4 A -1 , 
a considerable fraction of this Fermi surface is observed. 

Nevertheless, a more detailed analysis between the ex- 
perimental and theoretical data is hampered by the de- 
pendence of the ARPES results on the method of anal- 
ysis. Use of the maximal intensity at the Fermi level 
leads to Fermi surface contours, which compared to the 
theoretical predictions are shifted towards the occupied 
states. As a consequence, the corresponding electron-like 
Fermi surface is smaller, the hole-like Fermi surface big- 
ger than the LDA contours. In general, the deviations 
are small but the error depends on the absolute size of 
the Fermi surface and may reach —45 % in a few direc- 
tions. In contrast, analysis of the data with the help of 
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FIG. 16: (Color online) Cut through the Fermi body in the a*c*-plane of M0O2 as indicated in Fig. 1151 experimental Fermi 
surface contours from ARPES measurements on the (100) and (101) surfaces and theoretical Fermi surface contours from LDA 
calculations. 



the gradient of the photoemission intensity delivers con- 
tours, which agree much better with theory. However, 
again there exist a few single directions, where the devi- 
ations from theory may go up to —30%. Nevertheless, 
as the comparison of different methods of ARPES data 
analysis revealed, the gradient method delivers the best 
resultsi 52 ! 53 For this reason, only the results of the in- 
tensity gradient method will be included in the following 
discussions. 

Finally, in order to allow for a comparison with the 
de Haas-van Alphen data to be presented in the subse- 
quent section, the ARPES surface area of the electron- 
like Fermi surface in the b*c* -plane has also been de- 



termined and compared to the LDA results. Under the 
assumption of an ellipsoidal form the area can be calcu- 
lated using the Fermi vectors in the directions T-Y (see 
Fig. [TBI) an d r-Z (see Fig. [12]). As a result, the deviation 
of the theoretical surface area from the experimental one 
amounts to — 7% (see Tab. HIl below). 

D. De Haas- van Alphen (dHvA) oscillations 

The de Haas-van Alphen effect provides a well- 
established alternative to study the Fermi surface of 
metallic systems. This quantum mechanical effect 
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manifests itself as an oscillation of the magnetiza- 
tion/susceptibility as a function of an external magnetic 
field. The basic ideas of magnetic oscillations in metals 
were first described by Onsager^ and later on detailed 
by Shoenbergi^ 

The periodicity of the oscillations as function of the 
inverse of the magnetic field 1/|B| can be associated with 
a frequency F, which usually is given in units of Tcsla [T] . 
This frequency F is related to the extremal Fermi surface 
cross section Ap eTm i via the Onsager relation 



F[T] — A Fermi 

Zire 



(4) 



Variation of the direction of the magnetic field provides 
a good guess of the Fermi surface geometry. However, 
due to the ellipsoidal approximation, which is usually 
applied to derive the Fermi surface cut from the cross 
section ^Fcrmb the shape of the Fermi surface may not 
be fully resolved. Hence, for more complex Fermi surfaces 
the evaluation of the data must be complemented, e.g., 
by LDA calculations. In passing we point to the well- 
known fact that, in contrast to the ARPES experiments 
discussed above, de Haas-van Alphen measurements al- 
low to determine only the extremal cross sections of the 
Fermi body but do not have access to the complete band 
structure E(k). 

In a more quantitative treatment performed by Lifshitz 
and Kosevich^ and Shoenberg^ the oscillating part of 
the magnetisation parallel to B is given by the Lifshitz- 
Koscvich formula 
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4r4, Rt, Rd, and R$ are reduction 



where B = |B|, A" = 

coefficients, and the sum includes all harmonics p. 

While a closer investigation of the measured signal re- 
duction delivers more details of the electronic structure, 
in this paper only the temperature dependent factor Rt 
is of further interest. A finite temperature T causes a 
broadening of the Fermi energy, which is reflected by the 
reduction coefficient 



Rt = 



X 



sinh X 



X = 



2TT 2 pk B T 
fkjJc. 



(6) 



Finally, via temperature dependent measurements the ef- 
fective mass Tfi c of charge carriers included in u) r = — 
can be determined. 

For the dHvA measurements a 15 T magnet system in 
combination with a 3 He/ 4 He dilution refrigerator and a 
rotator of Oxford Instruments was used. In our setup the 
field modulation technique was used. Thus the dHvA 
magnetisation is a function of Bq + bo cos(io mo dt) and 
was detected via a mutual inductance bridge (astatic coil 
pair). The measurements were performed at B between 
10 and 15 T with a field variation rate of 15mT/min. 
The amplitude of the excitation field &o was a few mT 



with a modulation frequency of u> mo d — 19 Hz. With 
these values, heating of the specimen by eddy currents 
can nearly be neglected and the temperature was between 
15 and 25 mK. 



E. Results (dHvA) 

In order to determine the Fermi body of M0O2, the an- 
gular dependence of the dHvA oscillations was measured. 
With b* as the axis of rotation (which is perpendicular 
to the sheet plane of the following figures) the angle of 
the magnetic field was varied in the a* c* -plane in steps 
of 15°. Due to point symmetry the measurements ran 
only from 345° down to 150° (relative to the a-axis). To 
separate the oscillating part of the signal from the back- 
ground, a third order polynomial was subtracted from 
the 2048 points equidistant in 1/B (sec inset Fig. IT7| . A 
subsequent Fast Fourier Transformation yielded the fre- 
quency spectrum in units of Tesla (see Fig. [TT)) . After 



Oscillating part of the magnetisation 



III' (111 



j jfl'il'tfl W\ 



) I ,„l 




4000 6000 8000 10000 12000 

Frequency [T] 

FIG. 17: (Color online) dHvA frequencies of extremal cross 
sections of the Fermi body in one direction. All extremal 
cross sections are captured simultaneously (1531 T, 3943 T, 
4821 T), higher harmonics are given in brackets (3059 T, 
7885 T). The inset shows the underlying oscillating part of 
the magnetisation as function of 1/B after subtraction of a 
third order polynomial. 

elimination of higher harmonics and occasional isolated 
points the angular dependence of the frequencies - and 
the corresponding mirror points - was plotted in polar 
coordinates. Each point in the plot represents an area 
of extremal cross section of the Fermi body. 1 T corre- 
sponds to an area in k-space of 9.55T0 -5 A -2 . It must be 
pointed out that all extremal cross sections perpendicu- 
lar to B were captured simultaneously. This complicated 
the classification of different Fermi surface sheets and of 
possibly existing multiple extremal cross sections of one 
Fermi surface sheet. 

In the a*c*-planc of M0O2, three different Fermi sur- 
face sheets have been detected, which will be treated in 
detail in the following three subsections. 



1. dHvA measurement of sheet 1 (Y-hole) 
In Fig. [THJ the angular dependence of the dHvA fre- 




FIG. 18: (Color online) dHvA findings and LDA results of the 
hole-like Fermi surface around the Y-point (sheet 1) in the 
a*c*-plane. The size of the experimental points corresponds 
to the intensity of the signal. 

quencies of one Fermi surface sheet in the a* c* -plane is 
plotted. The correlation of the plotted points is well de- 
fined by the marginal change of frequency and the rela- 
tive high intensity. The latter shows a maximum between 
120° and 180° caused by a minimal curvature of the area 
(My oc -y=pO- The absence of appreciable intensity be- 
tween 30° and 75° is an indication of deviations of the 
rotational symmetry referred to the a-direction. The de- 
viations can be observed as a small bulge at an angle 
of 105°. The measured frequency amounts to 1126 T, 
whereas the frequency at the corresponding angle (with 
respect to the rotational symmetry) of 255° amounts 
only to 1076 T. This variation of the magnitude of the 
bulge in different directions (regarding the whole three- 
dimensional Fermi surface) can explain differences of the 
signal intensities. The extension of the dHvA Fermi sur- 
face in the c* -direction (T-Y) amounts to 1084 T, in di- 
rection perpendicular to c* to 1793 T. 

Shape and size of the experimental Fermi surface cross 
section allow the assignment to one of the LDA Fermi 
surfaces (see Fig. Uh)- The hole- like Fermi surface struc- 
ture centered around the Y-point (sheet 1) gives a very 
similar shape with a bulge at 100°. The size is likewise 
in good agreement with the results of the calculations. A 
quantitative comparison can be given for the c*-direction 
and the direction perpendicular to c*. The discrepancy 
of dHvA extremal surface areas to LDA calculations is 
— 11% in the c* -direction and —8% perpendicular to c*. 
Hence, the experimental dHvA Fermi surfaces are smaller 
than those theoretically predicted. 
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2. dHvA measurement of sheet 2 (Y-hole) 
In Fig. [19J the angular dependence of the dHvA fre- 
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FIG. 19: (Color online) Same as Fig. [18] but for sheet 2 (Y- 
hole). 

quencies of a second Fermi surface sheet in the a* c* -plane 
is shown. The shape of the dHvA Fermi surface resem- 
bles the hole-like LDA Fermi surface around the Y-point, 
called sheet 2 (sec Fig. Wfa)- However, the contraction of 
the experimental dHvA Fermi surface in direction of c* 
seems to be smaller than theoretically predicted. Here 
the deficiencies of the ellipsoidal approximation appear. 
Perpendicular to the a*c*-planc, in the b*c*-plane, the 
Fermi surface sheet 2 can not be described by an ellipsoid. 
In this plane, the LDA calculations yield a star shaped 
Fermi surface with an area larger than that of an ap- 
proximated ellipsoid (see Fig. |7fc). The deviations cause 
less contraction of the measured dHvA Fermi surface in 
c*-direction as predicted by theory. The high intensity 
of the signal at 90° indicates the continuation of the star 
shape of the b*c*-plane over the whole side. 

The measured frequencies amount to 2577 T and 
3205 T, parallel and perpendicular to the c*-axis, respec- 
tively. These results are —3% and —4% smaller than 
those predicted by the LDA calculations. 

3. dHvA measurement of sheet 3 (Electrons) 

Fig. EOl shows the angular dependence of the dHvA 
findings of a third Fermi surface sheet in the a*c*-planc. 
The shape is similar to the electron-like, peanut-shaped 
LDA Fermi surface sheet 3 (see Fig. UK)- Although the 
approximate rotational symmetry with respect to the a- 
axis, which was obtained by the LDA calculations (see 
Fig. Efc), simplifies the identification of the belonging 
points, there are large changes in frequency combined 
with weak intensity. The frequency of the peanut waist 
in direction of c* is 1166 T. The discrepancy to LDA cal- 
culations amounts to —14%. In the a-direction a fre- 
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FIG. 20: (Color online) Same as Fig. \TE\ but for sheet 3 
(electron pocket). 



quency of 5370 T is found. Here the deviation from the- 
ory amounts to —10%. Again, the experimental Fermi 
surface is smaller than the theoretical one. 

The missing signal at 315° and 330° can be explained 
by the open orbits caused by the neck of the peanut in 
the direction T-B. 



The quantitative comparison along selected high sym- 
metry directions shows that there are slight discrepan- 
cies between the dHvA and the LDA results. All dHvA 
Fermi surface areas are smaller than calculated by LDA. 
The discrepancies lie between —3% and —14%. The 
peanut shaped, electron-like Fermi surface sheet 3 shows 
the largest deviation to the LDA calculations. 

Due to the underlying ellipsoidal approximation, the 
shape of the dHvA Fermi surface sheets have to be 
treated with care. For some directions the real Fermi sur- 
face topology clearly differs from the assumed ellipsoidal 
shape. Especially, for sheet 2 (|| and _L to c*) and for 
sheet 3 (_L to c*) the LDA results indicate non- negligible 
deviations. 

A comparison of our dHvA results to previously 
published magneto-transport data reveals a very good 
accordancei 22 ' 23 i 29 ' 30 i 31 In particular, the results for the 
Fermi surface sheet 1 agree nearly perfectly, only differ- 
ences of less then 1.5% were found. In contrast, only 
little has been published for sheets 2 and 3. Yet, our 
results for these Fermi surface sheets, while being much 
more comprehensive, agree very well with the existing 
data, the differences being less than 6.5%. 



5. Effective masses 



4- Fermi surface in the a* c* -plane: dHvA vs. LDA 

In the following, the results of the dHvA measure- 
ments and the LDA calculations in the a*c*-plane will be 
compared. In Tab. [Til the dHvA surface areas and the 





||c* (90°) 


i-c* (0°) 


Sheet 1 


dHvA 


1084 T 


1793 T 




LDA 


1222 T 


1939 T 




(dHvA-LDA)/LDA 


-11% 


-8% 


Sheet 2 


dHvA 


2577 T 


3205 T 




LDA 


2650 T 


3333 T 




(dHvA-LDA) /LDA 


-3% 


-4% 


Sheet 3 


dHvA 


1166T 


5370 T 




LDA 


1358 T 


5968 T 




(dHvA-LDA) /LDA 


-14% 


-10% 




(ARPES-LDA) /LDA 


-7% 





TABLE II: Areas of Fermi surface cross sections as resulting 
from the de Haas-van Alphen measurements and the calcu- 
lations, respectively. The respective last lines for each sheet 
give the deviations of the experimental results from the cal- 
culated ones. For sheet 3, the corresponding deviation for the 
ARPES measurements has been added. 

LDA results parallel and perpendicular to c* are sum- 
marized. All three theoretically predicted Fermi surface 
sheets have been confirmed by the dHvA measurements. 
Qualitatively, the shapes of the experimental and theo- 
retical Fermi surface contours agree well. 



Temperature dependent measurements of the magnetic 
oscillations allow the determination of effective masses of 
the charge carriers via Eq. © . An increase in tempera- 
ture causes a reduction of the dHvA signal. For a single 
orientation of the magnetic field the effective masses of 
the three observed bands have been analyzed. The an- 
gle between B and a amounts to 75°. The results are 
displayed in Fig. 1211 The temperature has been varied 




o.o I rT7~TT i i i i i TT7~7-T-; 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Temperature [K] 



FIG. 21: (Color online) Temperature dependent measurement 
of the oscillation amplitudes for the three detected Fermi sur- 
face sheets and the corresponding effective masses using Eq. 
(|6|l as a fit. The angle between B and a amounts to 75°. 

between 0.1K and 3K. The oscillation amplitudes (nor- 
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malized to one at T = 0.1 K) as a function of the temper- 
ature are plotted. A fit of the data using Eq. © yields 
the effective masses m c . 

The effective masses of the three Fermi surface sheets 
differ clearly. Sheet 1 shows a slightly decreased ef- 
fective mass of 0.85m . However, sheets 2 and 3 ex- 
hibit increased masses of 1.87mo and 2.85mo, respec- 
tively. These values are in good agreement with existing 
Shubnikov-dc Haas measurements^ It has to be noted 
that the described geometry is not the one discussed in 
conjunction with the ARPES measurements. 



IV. DISCUSSION 

In general, we found excellent agreement between 
the experimental and theoretical Fermi surface contours. 
Both the ARPES and the de Haas-van Alphen results 
confirm the LDA calculations. In particular, all three 
Fermi surfaces sheets predicted by theory have been ver- 
ified experimentally, with respect to both the shape and 
the dimensions. 

Of course, while from the theoretical side all Fermi sur- 
face sheets can be calculated in one go, the measurements 
are restricted by the respective experimental setups and 
limitations. For instance, only certain cuts through the 
Fermi body are usually accessible. In addition, in the 
ARPES measurements matrix element effects may also 
play a role, which might undermine the experimental de- 
tection of portions of the Fermi surface. This may explain 
the fact that the calculated Fermi surface contours of all 
three sheets can only be partially observed by ARPES. 

While the overall agreement of the results of all three 
methods is very good, there still remain small differences 
on a quantitative level. These deviations have been in- 
cluded in Tabs. HI and HT1 To be specific, the deviations 
of the ARPES Fermi vectors from the calculated ones 
amount to several percent but may go up even to 30 % 
for some experimental points. Yet, both negative and 
positive deviations occur for both the electron- and hole- 
like surfaces, precluding systematic errors. This is dif- 
ferent for the de Haas- van Alphen measurements. All 
Fermi surface contours resulting from these experiments, 
including those given in the literature, are systematically 
smaller than the calculated onesi 22 i 23 i 29 ' 30 i 31 To be spe- 
cific, the electron-like Fermi surface sheet 3 is shifted up 
to —14% in direction of the occupied states, whereas 
the hole-like Fermi surface sheets 1 and 2 experience a 
shift of up to —10% in the direction of the unoccupied 
states. However, we point out that in general the differ- 
ences between the experimental and theoretical data are 
of a similar size as the differences between the different 
experimental techniques. 

Nevertheless, all these differences between the experi- 
mental and theoretical results, namely, both the unsys- 
tematic deviations of the ARPES data and the systemati- 
cally smaller dHvA Fermi surface contours, still conserve, 
within the above-mentioned limitations of the experimen- 



tal setups, the volume of the Fermi body, which, by con- 
struction, is correctly accounted for by the calculations. 57 

For the same reason, non-stoichiomctry of the samples 
can most likely be ruled out as a source of the observed 
quantitative differences. While the phase diagram for 
Mo-0 compounds excludes a hole-doping of more than 
1.3 a small defect concentration of 0.12% has been 
detected by Hall effect measurements.— In general, hole- 
doping causes a lowering of the Fermi level. According to 
the band structure shown in Fig. [6j this is equivalent to a 
shift of the Fermi surface contours towards the occupied 
states. However, a shift of up to —30 % would necessitate 
a hole-doping of up to 10%. Hence, non-stoichiometry 
will play only a minor role. 

Of course, one might argue, that electronic correla- 
tions may affect the quality of the calculated results. 
Indeed, our previous calculations infer the existence of 
weak electron correlations in M0O2' 11 ' 20 According to the 
present measurements, effective masses are only slightly 
enhanced. Yet, an enhancement may explain the reduc- 
tion of the band width as it has been observed in the 
k|| measurements on the MoO2(100) surface. However, 
magneto-transport measurements performed in this di- 
rection do not indicate an enhanced effective mass i 23 ' 24 i 29 
Furthermore, as has been demonstrated above part of the 
experimentally observed band width reduction could be 
accounted for by the inclusion of the full single-particle 
potential in the present calculations. 

Going into more detail, the deviations between the ex- 
perimental and theoretical data suggest that, as com- 
pared to the experiments, the calculated Fermi body is 
slightly shortened parallel to the a-axis and dilated per- 
pendicular to this axis. In the directions perpendicular 
to a this presumption is confirmed by the experimen- 
tal data. In particular, the cuts through the electron- 
like ARPES and dHvA Fermi body are smaller than the 
calculated results. The assumed shortening parallel to 
a is on somewhat weaker grounds. Only the waist of 
the "peanut" could be measured by ARPES, but no ex- 
perimental points exist at the top and bottom of the 
"peanut". In contrast, the dHvA measurements in the 
direction parallel to a might be affected by the under- 
lying ellipsoidal approximation. Yet, our assumption is 
well supported by the overall conservation of the volume 
of the enclosed Fermi surface as discussed above. Nev- 
ertheless, in passing we mention that these deviation are 
of the order of a few percent only. 



V. SUMMARY 

A detailed and comprehensive study of the electronic 
properties of M0O2 has been presented. This ma- 
terial has been analyzed using angle-resolved photo- 
electron spectroscopy (ARPES), de Haas- van Alphen 
(dHvA) measurements, and electronic structure calcula- 
tions based on a new full-potential implementation of the 
augmented spherical wave (ASW) method. The Fermi 
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surfaces determined by both kinds of experiments are in 
very good agreement with the theoretical predictions. In 
particular, three different Fermi surface sheets are cor- 
rectly identified by all three approaches. These sheets 
include an electron-like peanut-shaped Fermi surface cen- 
tered around the T-point, which has a volume of 11% of 
the Brillouin zone, and two smaller hole-like Fermi sur- 
faces both centered around the Y-point. 

As far as they can be detected within the limitations 
of the experimental geometries, slight differences between 
the epxerimental and theoretical results concern mainly 
the shape of the three Fermi bodies. To be specific, from 
the calculations, the peanut-shaped Fermi body appears 
to be slightly shortened and elongated parallel and per- 
pendicular to the a-axis. 

With the help of the new calculations, controversies 
concerning two additional Fermi surface sheets centered 
around the Z- and B-point could eventually be resolved. 
In particular, the occurence of these sheets could be at- 
tributed to the atomic-sphere approximation, which was 
used in our previous work. In contrast, in the new full- 
potential ASW calculations, these additional sheets do 
no longer appear. Furthermore, the full-potential treat- 
ment led to a much improved agreement between the 
ARPES results and the calculations as concerns the oc- 
cupied band width at the T-point. 

In general, the very good agreement between the re- 
sults obtained by the ARPES and dHvA measurements 
for metallic M0O2 has important consequences for the in- 
vestigation of those neighbouring transition-metal diox- 
ides, which display a metal-insulator transition. While 



the Fermi surfaces of the metallic phases of these mate- 
rials are not accessible to a low-temperature technique 
as dHvA experiments, the good agreement found in 
the present study implies that ARPES measurements 
will be useful to study the Fermi surfaces of the high- 
temperature metallic phases. 

Furthermore, the good agreement between the differ- 
ent experimental techniques and the electronic structure 
calculations, with the deviations between measured and 
calculated data being of the same order as those between 
the different experimental techniques, support our orig- 
inal proposal that these calculations are well suited to 
study the electronic properties of the early transition- 
metal dioxides. In particular, our results are very en- 
couraging concerning the validity of our calculations for 
the d 1 members VO2 and NbC>2 and support the inter- 
pretation of the metal-insulator transitions of the latter 
two compounds as embedded Pcicrls-likc instabilities. 
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